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Abstract 

We compute the b quark mass from dynamical lattice QCD with clover 
quarks. The calculation is done at a fixed lattice spacing with sea quark 
masses as low as half the strange quark mass. Our final result is m^m^) 
= 4.25(2) (11) GeV, where the first error is statistical and the last error is 
the systematic uncertainty. 

1 Introduction 

The mass of the bottom quark is a fundamental parameter of the standard 
model [1]. To extract the b mass from experiment, QCD corrections must be 
computed reliably. The best way to do this is use lattice QCD. The different 
methods of computing the mass of the bottom quark have recently been reviewed 
by El-Khadra and Luke [2]. 

The particle data table quotes the mass of the b quark in the MS scheme 
at the b mass to lie between 4.0 and 4.5 GeV [3]. It is particularly important 
to reduce the error on the b quark mass, because it is the cause of the largest 
uncertainty on the determination of V u b from the total inclusive B meson decay 
b — + ulv. El-Khadra and Luke [2] note that a 100 MeV error on nib corresponds 
to a 6 % error on the V u b CKM matrix element (currently only known to 19 % 
accuracy [3]). 

In this paper we use unquenched lattice QCD with static-light mesons to 
extract nib- As we discuss in section 4, the error due to the use of static (leading 
order HQET) approximation is only of order 30 MeV [4], hence the static limit 
has an important role for the phenomenology of determining V u b- 

2 Details of lattice calculations 

We used non-perturbatively improved clover fermions in both the sea and va- 
lence quarks. The Wilson gauge action is used for the gluons. The full details of 
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Name 


No. 


romps 


^sea 


™val 


Volume 


a£ 


^static GeV 


DF1 


20 


1.92(4) 


0.1395 


0.1395 


12 3 24 


0.87(1) 


0.59(5) 


DF2 


78 


1.94(3) 


0.1395 


0.1395 


16 3 24 


0.842(5) 


0.55(2) 


DF3 


60 


1.93(3) 


0.1350 


0.1350 


16 3 32 


0.772t£ 


0.69(7) 


DF4 


60 


1.48(3) 


0.1355 


0.1355 


16 3 32 


0.739±1 


0.66(5) 


DF5 


60 


1.82(3) 


0.1355 


0.1350 


16 3 32 


0.748±| 


0.68(5) 


DF6 


55 


1.06(3) 


0.1358 


0.1358 


16 3 32 


0.707ti2 


0.64(7) 



Table 1: Lattice binding energy (a£) and physical binding energy {h. static) for 
each data set. 



the actions and details of the unquenched calculation are described in [5, 6, 7]. 

We use static quarks for the heavy mesons. The lattice binding energy is 
extracted from a two point correlator using variational smearing techniques. 
The local two point function is 

C(t) = ^(0 | $ B (x,t)& B (x,0) | 0) (1) 

X 

= Z 2 cxp(-a£t) (2) 

where $b is the interpolating operator for static-light mesons. We have already 
published [8] an extensive analysis of the spectrum of static-light mesons. Our 
previous paper [8] also describes the all-to-all propagators used to improve the 
statistical accuracy and the fuzzing methods used. 

In table 1 we present our results for the lattice binding energy All the 
data sets used (i = 5.2. Data sets DF1 and DF2 used a clover coefficient of 
1.76, while all the others used the non-perturbative value of 2.0171. The results 
for the data sets: DF1, DF2 have already been published [8]. The ensemble 
size for data sets DF4 and DF5 have been trebled over the results previously 
published [8]. The data from ensembles DF5 and DF6 are new. 

3 Extracting the quark mass 

We evaluate the mass of a pseudoscalar heavy-light meson from lattice QCD 
with static heavy quark and compare with the experimental mass value. This 
gives information about the 6-quark mass. The strange quark mass is accessible 
in lattice evaluations, so to minimise extrapolation, we use the B s meson for this 
comparison. We still need to extrapolate the sea quark mass to the experimental 
value and we discuss this later. 

In this section we will describe the central values for our calculation. We 
discuss systematic uncertainties in section 4. The quantity £ , from the lattice 
calculation, contains an unphysical -j divergence (5m) that must be subtracted 
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off to obtain the physical binding energy (Astatic)- 

Astatic = £ -5m (3) 

The pole quark mass is determined from 

ml ole = M Bs - A statlc (4) 

The physical value [3] of the meson mass Mb s (5.369 GeV) is used. 

In the static theory 5m has been calculated to two loops by Martinelli and 
Sachrajda [9]. 

a5m = 2.1173a s (m&) + {(3.707- 0.225n/) log(m£a) 

- 1.306 -n/(0. 104 + 0Ac sw - 0A03c 2 sw )}a s (rnE) 2 (5) 

where n/ is the number of sea quarks and csw is the coefficient of the clover 
term. We discuss estimates of the next order to a5m in section 4. 

The pole mass (see Kronfeld for a review [10]) is converted to MS using 
continuum perturbation theory [11]. 

mp(M) - Z pm (v)ml ole + 0(1 /m b ) (6) 

where 

Z pm (^ = m b ) = 1 - (11.66 - 1.04n/)( ) z (7) 

The perturbative correction between the pole mass and the MS mass is known 
to 0(a 3 ) [12, 13]. The perturbative series connecting the pole mass with the 
MS mass is badly behaved due to renormalons (see [14] for a review). The 
lattice matching is only done to 0(a 2 ), hence we convert the pole mass to MS 
at the same order, using a consistent coupling, so the differences in the series 
are physical. 

We use the values of the coupling using the values of Aqcd from the joint 
UKQCD and QCDSF paper [15]. The four loop expression for a s is used [16] 
to determine the coupling from Aqcd- We consistently use n/ = 2 in all the 
perturbative expressions. For n sea = 0.1355 (0.1350) we use Aqcd = 0.178 ( 
0.173) MeV [15]. We use the same value of Aqcd for the two data sets DF1 
and DF2, where Aqcd has not been computed. 

In [8] we estimated the mass of the strange quark using the pseudoscalar 
made out of strange quarks [17]. This provided r mps = 1.84. This value is 
close to r mps — 1.82 jl^ for DF5 data set [6] hence we use the binding energy 
from that data set as the value at strange. This is a partially quenched analysis. 
The data sets DF1 and DF2 also have a sea quark mass close to the strange 
quark mass [17]. 
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To determine the lattice spacing we use the measured value for ro/a from 
the potential with the 'physical' value of ro as 0.525(25) fm [18]. We discuss 
in more detail the systematic error from the choice of ro in section 4. Hence 
our best estimate of mb(mb) — 4.25(2) GeV from the DF5 data set, where the 
errors are statistical only. 



4 Computing the systematic uncertainties 

Gimenez et al. [4] discuss the systematic error from the neglect of the l/m& 
terms in the static limit. Heavy quark effective field theory parametrises the 
heavy mass corrections to the mass of a heavy- light meson Mb ■ 

M B = ™ fc + A s ^ c -A--gL (8) 

where A sta tic is the static binding energy, Ai is the matrix element due to the 
insertion of the kinetic energy and A2 is the matrix element due to the insertion 
of the chromomagnetic operator. The value of A2 ~ 0.12 GeV 2 can be obtained 
from the experimental mass splitting between the B* and B mesons. The value 
of Ai is much harder to estimate. Gimenez et al. [4] use a range of Ai from -0.5 
to 0.0 in GeV 2 . This includes the determination from quenched lattice QCD of 
Ai = -(0.45 ± 0.12)GcV 2 by Kronfcld and Simonc [19]. JLQCD have recently 
tried to compute Ai using NRQCD [20]. As suggested by Gimenez et al. [4], 
using a symmetric error of 30 MeV to parameterise the neglected 1 / mj terms 
seems reasonable to us. 

The specification of the strange quark mass, described in section 3, essen- 
tially relies on the physical K mass. Since the <f> meson has a very narrow 
width, it may also be used to specify the strange quark mass. In their study 
using essentially the same lattice parameters, JLQCD [21] see approximately a 
10% difference between using the <j> and the K to set the strange quark mass. 
Motivated by JLQCD's result, we use a symmetric error of 5% as an estimate 
of the additional uncertainty in our estimate of the strange quark mass. This 
induces an error of 60 MeV in mb(mb) for the DF5 data set. 

The data sets DF1, DF2 were generated using a different value of csw to 
that used to generate data set DF5, hence they can not be used to estimate the 
size of the lattice spacing effects. The comparison of the results between data 
set DF1 and DF2 can in principle be used to estimate finite size effects. As 
the physical size of a size of the lattice changes from 1.83 fm (DF1) to 2.44 fm 
(DF2), mb(rrib) changes from 4.33(2) to 4.369(7) GeV. Hence, a simple estimate 
of the finite size effects in date set DF5 (size of box 1.77fm) is -39 MeV. In 
quenched QCD Duncan et al. [22] found no finite size effects in the binding 
energy for lattice lengths: 1.3, 1.8 and 2.2 fm. Although, finite size effects in 
quenched and unquenched QCD can be very different, we think it more likely 
that the differences in two data sets is due to a statistical fluctuation on the 
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smaller lattice. This is supported by the fact that UKQCD saw no finite size 
effects in the light hadron spectrum between DF1 and DF2 [5]. 

The choice of coupling is a systematic error. Gimenez et al. [4] use J^qcjj 
= 300 MeV as the central value. We use the result for Aq C d determined from 
the DF4: and DF5 data sets [15]. We don't feel that it is appropriate to use 
the values of Aqcd from experiment (as done by Gimenez et al. [4]) even for 
deriving a systematic error. The agreement between Aqcd from lattice QCD 
calculations and experiment is not good for calculations that use clover fermions 
for the sea quarks [21, 23, 15]. We assume that the discrepancy will be reduced 
as calculations are done with lighter sea quark masses and finer lattice spacings. 

To estimate the effect of the chiral extrapolation of the sea quark masses, 
we extrapolated A sta tic from data sets DFA and DF5 linearly in (romps) 2 to 
r o m PS — 1-93 (the same as for data set DF3). The extrapolated result for 
Astatic at r m PS = 1.93 at n sea = 0.1355 was consistent with A static on data 
set DF3 (Ksea — 0.1350). We see no evidence for the dependence of A sta tic on 
the sea quark mass. Gimenez et al. [4] see a slight increase in the lattice binding 
energy with decreasing quark mass. There are potentially non-analytic m PS 
terms in the mass dependence of the binding energy [24] . 

To estimate the systematic errors on the perturbative matching we did a 
number of things. Following Gimenez et al. [4] we compared taking the product 
of the two perturbative factors in equation 6 against expanding the perturbative 
expressions and only keeping 0(a 2 s ) terms. This increases the mass for data set 
DF5 by 25 MeV. 

In quenched QCD the next order correction to equation 5 has been computed 
numerically by two groups using different techniques [25, 26]. As the two groups 
obtained essentially the same result we will use the result of Di Renzo and 
Scorzato [25]. In quenched QCD the next order correction to equation 5 is [4] 

aSm {3) = (A^+6.48945(-3.57877+log(mFa))(3.29596+log(Wa)))a ;s (7n^) 3 (9) 

where Xi is the number from the numerical calculation. Di Renzo and Scorzato [25] 
obtain X2 = 86.2(0.6)(1.0) for quenched QCD. Di Renzo and Scorzato have com- 
puted X 2 for ri/ = 2 Wilson fermions [27]. The new result also involves a lattice 
calculation of the MS coupling for Wilson fermions, so it is not obvious how to 
incorporate the new result into this analysis. 

The next order to the connection between the pole mass and MS mass is 
known in the continuum (equation 7). 

Zfl = -(157.116 - 23.877971/ + 0.6527n2)(^^) 3 (10) 

Because equation 9 is known for quenched QCD, we do not use it for our central 
result. It does seem appropriate to use equation 9 and equation 10 to estimate 
the systematic errors due to the neglect of higher order terms. We set n/ = 
in equation 10. Adding in the next order term reduces the mass of the bottom 
quark by 12 MeV for data set DF5. 
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Because of limitations in computer time, the unquenched calculations are at 
fixed lattice spacing (so the continuum limit has not been taken) and fairly heavy 
sea quark masses are used. This means that different lattice quantities produce 
slightly different values of the lattice spacing. In [18] the values of ro from various 
calculations that used clover fermions are collected together. The results were 
in the range r = 0.5 to 0.55 fm. This motivates our choice of ro = 0.525(25) 
fm. Using mass splittings in Upsilon on improved staggered configurations with 
measurements of the potential M1LC [28] quote r = 0.467 fm. This is 2.3a from 
our central value for ro. In the graph [29] of spread of variation of lattice spacings 
from different physical quantities, the most dramatic failures of the quenched 
approximation occur for P — S and 25 — IS mass splittings in Upsilon and the 
pion decay constant. We speculate that these quantities are more sensitive to 
the heavy quark potential at the origin that depends on n/ from the running of 
the coupling (see [30] for a discussion of this). The relatively strong dependence 
of the Upsilon mass splittings and pion decay constant on the sea quark mass 
and lattice spacing doesn't make them a good choice to set the scale for current 
unquenched calculations with clover fermions. Hence we feel that using the 
value from M1LC for ro (as advocated in [31]) will artificially inflate the error 
bars, so we stick to our original estimate of r = 0.525(25) fm. 

The perturbative analysis, used in the this section, assumes that the sea 
quark mass is zero. However, the sea quark masses used in current calculations 
with Wilson-like quarks arc not negligible. At n sea = 0.1355 and 0.1350, the 
vector definitions of the sea quark mass (in units of the lattice spacing) are 0.026 
and 0.044 respectively. The light quark mass dependence has been computed 
by Bali and Boyle [32]. 

For light quark masses below 0.1 (in lattice units) Bali and Boyle [32] provide 
a quadratic parameterisation of the light quark mass dependence of 5m. The 
expression in equation 5 gets modified to 

aSm = 2.1171a a (m&) + ((3.707 - 0.225n/) log(m^a) 

- 1.306 -n / (-0.199 + 0.516m 9 -0.421m2))a iS (mF) 2 (11) 

where we have specialised to csw = 1- Equation 5 is a combination of the two 
loop static self energy and a conversion from the bare coupling to the massless 
MS scheme [9]. As stressed by Martinelli and Sachrajda [9], it is important to 
use a consistent coupling in equation 6, so that the poorly behaved perturbative 
expansion of Z pm cancels with that of 6m. This makes it easier to use the 
massless quark MS scheme. The determination of Aqcd includes the effects of 
the masses of the light quarks [15]. The use of 11 changes the central value of 
mb(rrib) by 1 MeV for the data set DF5. 

For our final result we use the central value from data set DF5. The sys- 
tematic uncertainties have been discussed in this section. Hence our final result 
is 

m{m) = (4.25 ± 0.02 ± 0.03 ± 0.03 ± 0.08 ± 0.06)GeV (12) 
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Group 


comment 


mb(mb) GeV 


This work 
Collins [35] 
Gimenez et al. [4] 
Di Renzo and Scorzato [4] 

Bali and Pineda [36] 
Heitger and Sommer [34] 


unquenched 
unquenched 
unquenched 
unquenched 
quenched 
quenched 


4.25(2)(11) 
4.34(7)1? 
4.26(9) 
4.21 ±0.03 ±0.05 ±0.04 
4.19(6)(15) 
4.12(8) 



Table 2: Lattice QCD results for mb(rrib) in the MS scheme. The last error on 
the Bali and Pineda result is an estimate of unquenching. 



where the errors are (from left to right): statistical, perturbative, and neglect of 
l/mft terms, ambiguities in the choice of lattice spacing, and error in the choice 
of the mass of the strange quark. 
Gimenez et al. [4] obtain 

mE(mE) = (4.26 ± 0.03 ± 0.05 ± 0.07)GeV (13) 

from a simulation at fi = 5.6 and volume = 24 3 40 with two dynamical quark 
masses, from the T\L collaboration. The first error is due to statistics. The 
second error includes the neglect of the 1 /rrib terms and the ambiguity in the 
determination of the lattice spacing. The third error is due to the neglect of 
higher order corrections in the perturbative matching. Gimenez et al. [4] used a 
preliminary result for X 2 , that was relatively imprecise [33] hence their estimate 
of the higher order effects is looser than ours. Also we used Aqcd determined 
consistently from this data set, while Gimenez et al. [4] used continuum based 
estimates of the coupling. This analysis has recently been updated by Di Renzo 
and Scorzato [27] with the unquenched value of X 2 - Gimenez et al. [4] only used 
one quantity to estimate the lattice spacing, so their error from the ambiguity 
in the choice of lattice spacing is underestimated in the final result. However, 
the two effects compensate and the final error is probably representative. It 
is pleasing that our calculation with a different set of parameters is essentially 
consistent with that of Gimenez et al. [4]. 

5 Conclusions 

In table 2 we collect some recent results for the mass of the bottom quark from 
lattice QCD. Our result is consistent with the previous unquenched calculations. 
Unfortunately, we have not managed to reduce the size of the error bars. The 
largest error in the recent values for the mass of the bottom quark is due to 
the spread in different lattice spacings. Heitger and Sommer [34] noted that a 
change in r by 10% changed the value of m^rrib) by 150 MeV. 

The prospects for an improved estimate of the mass of the bottom quark 
from lattice QCD are quite good. The unquenched calculations with improved 
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staggered quarks produce consistent lattice spacings from many different quan- 
tities [29]. The numerical calculation of the third order contribution to Sm 
has been done for unquenched Wilson fermions [27]. Applying the technology 
of Heitger and Sommer [34] to unquenched calculations would allow a non- 
perturbative estimate of the mass of the bottom quark that is free from prob- 
lems with delicate cancellations of poorly converging perturbative expressions. 
The use of automated perturbative calculations may allow the computation of 
the bottom quark mass from NRQCD / FNAL type calculations with two loop 
accuracy [37, 38]. We are investigating the use of the static formulation in- 
troduced by the ALPHA collaboration [39], but the required perturbative (or 
non-perturbative) factors are not available yet. 

Some combination of the techniques and projects mentioned in the last para- 
graph should be able to produce a number of independent calculations, with 
different systematic errors, of the mass of the bottom quark from unquenched 
lattice QCD. 
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